Polynomial perturbations of hermitian 
linear functionals and difference equations 

M.J. Cantero, L. Moral, L. Velazquez * 
Departamcnto de Matematica Aplicada. Univcrsidad dc Zaragoza. Spain 



Abstract 

This paper is devoted to the study of general (Laurent) polynomial 
modifications of moment functionals on the unit circle, i.e., associated 
with hermitian Toeplitz matrices. We present a new approach which 
allows us to study polynomial modifications of arbitrary degree. 

The main objective is the characterization of the quasi-definiteness 
of the functionals involved in the problem in terms of a difference 
equation relating the corresponding Schur parameters. The results are 
presented in the general framework of (non necessarily quasi-definite) 
hermitian functionals, so that the maximum number of orthogonal 
polynomials is characterized by the number of consistent steps of an 
algorithm based on the referred recurrence for the Schur parameters. 

Some concrete applications to the study of orthogonal polynomi- 
als on the unit circle show the effectiveness of this new approach: an 
exhaustive and instructive analysis of the functionals coming from a 
general inverse polynomial perturbation of degree one for the Lebesgue 
measure; the classification of those pairs of orthogonal polynomials 
connected by a kind of linear relation with constant polynomial coef- 
ficients; and the determination of those orthogonal polynomials whose 
associated ones are related to a degree one polynomial modification of 
the original orthogonality functional. 
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1 Introduction 



The intense activity during the last decades around the theory of orthogonal 
polynomials on the unit circle has stimulated the study of perturbations of 
hermitian functionals. The possibility of considering modifications that do 
not preserve the hermitian character of the functional leads to left and right 
orthogonal polynomials (see [3]), thus most of the efforts have been concen- 
trated in the analysis of hermitian perturbations as a source of new families 
of standard orthogonal polynomials (see the recent monograph on orthogo- 
nal polynomials on the unit circle \^22\ [23] ) and the references therein) . 

This paper proposes a new method to study the hermitian modifications 
obtained when multiplying a hermitian functional by a Laurent polynomial 
of any degree, in short, the hermitian polynomial modifications. This kind 
of perturbation has been considered previously (see for instance [H HI [U El [71 
[g [1 [M El [El [II [13 [El [13 [ig [11 [M] ), but the usual approaches have the 
drawback of being formulated in terms of orthogonal polynomials, kernels 
and determinants, what makes difficult the practical application, specially 
for perturbations of high degree. 

On the contrary, our method is based on a recurrence for the Schur pa- 
rameters of the two functionals involved in the perturbation. This provides 
an algorithm to generate the Schur parameters of one of the functionals, 
starting from the Schur parameters of the other functional. Furthermore, 
this recurrence yields a characterization of the maximum number of orthog- 
onal polynomials for one of the functionals, given the number of orthogonal 
polynomials that the other functional has. That is, this approach permits 
us to study the relation between the quasi-definiteness of a functional and 
a polynomial modification of any degree with less computational effort than 
the methods already existing. 

We distinguish between three different but related problems, depending 
on the data at hand: 

• Basic problem: characterize when two functionals are related by a 
polynomial perturbation in terms of their Schur parameters. 

• Direct problem: characterize the quasi-definiteness of a polynomial 
modification from the Schur parameters of the original functional. 

• Inverse problem: characterize the quasi-definiteness of a functional 
from the Schur parameters of one of its polynomial modifications. 

Despite the symmetry between the direct and inverse problems, they 
have a quite different nature which makes much more interesting the last 
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one. The root of this difference is the fact that, given a functional and a 
Laurent polynomial, the corresponding polynomial modification is uniquely 
defined while there are infinitely many functionals whose modification is 
the given one. This leads to a rich structure in the set of solutions of the 
inverse problem which, as we will see, is related to another kind of interesting 
modifications: the addition of Dirac deltas and its derivatives. Hence, any 
information about inverse polynomial modifications can be translated as a 
result on perturbations by Dirac deltas. 

Furthermore, as an example will show, some special solutions of an in- 
verse problem can act as "attractors" for the asymptotics of the parameters 
of other solutions. Thus, the analysis of those special solutions provides 
information about the asymptotics of perturbations by Dirac deltas. 

The rich structure of the inverse problem has a double interest due to the 
fact that our approach, based on a recurrence for the Schur parameters, also 
yields interesting connections between the study of polynomial modifications 
and difference equations. Therefore, the asymptotics of the solutions of the 
inverse problem is closely related to the asymptotics of difference equations. 

The content of the paper is structured in the following way: the rest 
of the introduction summarizes the basic definitions and notations; Section 
[2] includes the main results about hermitian polynomial modifications, i.e., 
it is devoted to what we call basic problem; direct and inverse problems 
are discussed in Section [Sj including an exhaustive analysis of an explicit 
example of inverse problem; and Section U] shows other applications of the 
techniques developed in the paper, i.e., a complete classification of the pairs 
of orthogonal polynomials related by certain type of linear relations with 
constant polynomial coefficients, and the determination of the orthogonal 
polynomials whose associated ones come from a polynomial modification of 
degree one of the original orthogonality functional. 

Now we proceed with the conventions for the notation. 

In what follows T := {z e C : \z\ = 1} and D := {2 e C : |z| < 1} 
are called respectively the unit circle and the open unit disk on the complex 
plane. P := C[z] is the complex vector space of polynomials with complex 
coefficients, and the vector subspace of polynomials whose degree is not 
greater than n, while P_i := {0} is the trivial subspace. A := C[z, z~^ is 
the complex vector space of Laurent polynomials and, for m < n, we define 
the vector subspace Am.n '■= span{z"^, . . . , z"}. Given any / € A 

we define = f{z~^) and, if p G P^ \ Pn-i, P* denotes its reversed 

polynomial = z"'p^{z). Sometimes we use the notation = z'^p^,{z) 
for polynomials p G P„ whose degree can be smaller than n. Then we refer 
to the *n operator when it is advisable to avoid misunderstandings. 
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Any hermitian linear functional f on A {v[z "] = n = 0,1,...) 

defines a sesquilinear functional (• , •)^: A x A — > C by 

{f,g)v-=v[f*g], f,geA. 

The sequence {pn)n>o is a sequence of orthogonal polynomials with respect 
to the hermitian linear functional v if 

(i) p„GP„\P„_i, 

(ii) {pniPm)v — ^n^n,mi 7^ 0, 

and when such a sequence exists v is called a quasi-definite functional. If 
v[l] 7^ we can assure only the existence of a finite segment of orthogonal 
polynomials, i.e., a finite set {pk)^=Q of polynomials satisfying (i) and (ii). 
When V has a finite segment of orthogonal polynomials {pk)k=Q of length 
n + 1 we say that v is quasi-definite on P„. 

In the positive definite case (/„ > 0, n = 0, 1, . . .) there exists a posi- 
tive measure /i supported on T providing an integral representation for the 
functional 

v[f] = I f{z)dfi{z), /eA. 

Due to this reason a sequence {pn)n>o satisfying (i) and (ii) is called a 
sequence of orthogonal polynomials on the unit circle, even in the general 
quasi-definite case. If In = ±1 for all n, (p„) is called a sequence of orthonor- 
mal polynomials on the unit circle. We denote by (j5„)„>o the orthonormal 
polynomials with positive leading coefficients. 

In that follows {ipn)n>o denotes the sequence of monic orthogonal poly- 
nomials (MOP) with respect to a hermitian functional v. Two hermitian 
linear functionals vi,V2 have a common finite segment (V'j)j=o of MOP iff 
there exists A E M* such that vi[f] = \v2[f] for any / E A_„^„, although 
requiring this condition to hold only for any / E Pn is enough due to the 
hermiticity. In this case we say that vi and V2 are equivalent in P„ or, in 
a more symbolic way, vi = V2 in P„. If this holds for any n, we simply say 
that vi and V2 are equivalent and we write vi = V2- 

A sequence {'tpn) is a sequence of MOP on the unit circle iff it satisfies 
the recurrence relation (see [25l [HI [22] ) 

tpniz) = Z'il:n-i{z) +'il:n{^Wn-i{z)^ n = 1, 2 . . . , (1) 

with ipo{z) = 1 and 1^/^^(0)1 ^ 1 for n > 1. Applying the operator to the 
above recurrence we get the equivalent one 

rn{z)=ihMz^n-l{z)+rn-l{z), n = 1, 2 . . . . (2) 
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The values tpn{0) are called the Schur parameters or reflection coefficients 
of the hermitian linear functional v. 

A straightforward computation yields 

1-|^„(0)|2 = ^, n=l,2,..., 

where e„ := (V'n,V'n)«; = v[ipnZ~'^] relates p„ and tpn by p„ = 

When V is positive definite e„ = ||V'n|||2(^) > for n > 0, which means that 

\^n{0)\ < 1 for n > 1. 

2 Hermitian polynomial modifications 

We are interested in those (Laurent) polynomial modifications of hermitian 
functionals which preserve their hermitian character, in short, the hermitian 
polynomial modifications of hermitian functionals. If u is a linear functional 
on A and L G A the modified functional vL is defined by 

vL[f]:=v[Lf], feA. 

The modified functional vL is hermitian for every hermitian t; iff L* = L, 
which is equivalent to state that L = P + with P S P (see [2]). Such a 
polynomial P can be uniquely determined by L simply requiring P(0) G M, 
a convention that we will assume in what follows. We will refer to degP as 
the degree of the polynomial modification, which we will consider greater 
than or equal to one, and L will be called a hermitian Laurent polynomial 
of degree r. 

Another way to characterize a hermitian polynomial modification is 
through the polynomial A = z'^'^^^L of degree 2degP. The condition 
means that A is self reciprocal, i.e.. A* = A. Thus the hermi- 
tian polynomial modifications are related to the self-reciprocal polynomials 
of even degree. 

The set of roots of a self-reciprocal polynomial, counting the multiplicity, 
is invariant under the transformation Q 1/C- That is, their roots lie on 
the unit circle or appear in symmetric pairs C)l/C- Indeed, this property 
characterizes the self-reciprocal polynomials up to numerical factors. This 
implies that any self-reciprocal polynomial of even degree factorizes into a 
product of self-reciprocal polynomials of degree 2. As a consequence, an 
arbitrary hermitian polynomial modification is a composition of elementary 
ones of degree 1, i.e., if L = P + P* with degP = r, then L = L1L2 • • • Lr 
with Lk = + Pfc^, and degP^ = 1. 
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Sometimes we will deal with polynomials A G P„ whose degree is not 
necessarily n but such that A*" = A. In this case we will say that A is 
self-reciprocal in IPfi to avoid misunderstandings. Such a polynomial has 
the general form A(z) = z^B{z) where B is strictly self-reciprocal. Thus, a 
self-reciprocal polynomial in P„ is actually self-reciprocal iff it has no zeros 
at the origin. 

Given a hermitian functional v and a Laurent polynomial L = P + P^, 
our purpose is to obtain relations between the MOP and Schur parameters 
associated with the functionals v and vL. Multiplying L by a non null 
real factor gives rise to a hermitian functional which is equivalent to vL 
and, hence, with the same MOP and Schur parameters as vL. Therefore, 
concerning our aim, the Laurent polynomial L, as well as the polynomials 
P and A, are defined up to non null real factors. 

The following general result will be useful to achieve our objective. In 
what follows we denote by S'^" the orthogonal complement in IP-;^ of a sub- 
space 5 C Pn- 

Lemma 2.1 (see [26]). Let v be a hermitian functional such that the cor- 
responding n-th MOP ipn exists. Then, 55 = {z^il!nYk=Q ^ {'^^'^nYk^o ^ 
basis of (z^Fn-r-i)^"^'' for n > r > 1, and a generator system o/P^+r for 
r > n > 0. 

Sketch of the proof. If n > r > 1, the orthogonality of ipn assures that 53 C 
(z^P„_,._l)-^"+^ Besides, n E span«B iff = Ci:n + Dijj*^, C G P^^, D G P,._i. 
Furthermore, this decomposition is unique because gcd(V'n5V'n) = which 
proves the linear independence of ®. Then, the first result follows from the 
fact that tj«B = 2r + 1 = dim(z'^P„„r_i)^"+". 

Suppose now that r > n > 0. From the previous result we know that 

{^Hn}l=o U {z'^rnVkll is a basis of P2n. Heuce, {z'^VnlLo U {^'^^;}fc=d 
is a linear independent subset of Pn+r with n + r + 1 elements, thus it is a 
basis of Pn+r 5 which proves the second result. □ 

Our interest in the previous lemma is the following direct consequence. 

Corollary 2.2. Let v be a hermitian functional such that the corresponding 
n-th MOP ifjn exists. Then, every polynomial Q, € (z''P„_r-i)^"^'^ has a 
unique decomposition VL = C^l^n + L)ipn, C G Fr, D £ Pr-i, for n> r > 1, 
and every polynomial G Pn+r has infinitely many such decompositions for 
r > n > 0. 

Remark 2.3. It is worth it to remark the case n = r in the above corollary, 
which says that every polynomial $7 G P2r admits a unique decomposition 

n = Cipr + Dip*, C €Fr, D e Fr-l. 
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The next theorem is the starting point for our approach to the study of 
hermitian polynomial modifications of hermitian functionals. 

Theorem 2.4. Let u, v be hermitian functionals with finite segments of 
MOP {(Pj)J=o, [ipjTjto respectively, and let L = P + = z"'' A with P a 
polynomial of degree r. Then, the following statements are equivalent: 

(i) u = vL in P„. 

(a) There exist Cj G Pr, Dj € Pr-i with Cj{0) 7^ such that 

Aifj = Cjipj+r + Dj'il)*^^^, j = 0,...,n. (3) 

(Hi) There exist Cj G Pr, Dj G Pr-i with Cj{0) 7^ such that 

A^] = zD*^,+r + C*rj+r, D* = D*f~\ j = 0,...,n. (4) 

The polynomials Cj G P^, Dj G Pr-i satisfying (0) or Q are unique, 
deg Cj = r, Cj{0) G M and C*(0) = ^(0). 

Proof. The equivalence between (ii) and (iii) follows from the use of the *2r+j 
operator and the fact that A is a self-reciprocal polynomial of degree 2r. 
Also, assuming (ii) we get degCj = r because deg(DjVj+r) < deg(^v'j) = 
2r + j, and the equality {ipj,Lpj)u = u[ipjZ~^= Cj{0)ej+r implies Cj(0) G M. 
On the other hand, evaluating ([4]) at z = we find that C*(0) = ^(0). It 
only remains to prove the equivalence between (i) and (ii) and the uniqueness 
of decomposition ([3]). 

Suppose (i), i.e., u[f] = XvL[f], A G M*, for any / G A_„^„. The 
orthogonality of {ipj)'j^Q with respect to u gives 

= u[ipjZ~''] = Xv[AipjZ~'^''~^'''^], r<k + r<j + r-l, 

which means that Aipj G {z'''¥j-i)'^^''^^ with respect to v. Using Corollary 
12.21 we get ([3]) and the uniqueness of the polynomials Cj , Dj . 

On the other hand, if {y^j)^=Q, {ipj)^^Q satisfy ([3]), the orthogonality of 
i'^j)^=o with respect to v yields 

vL[^jZ-'']= v[AipjZ-^^+'^]= v[iC,^P,+r + I?jV;+.)^"^'+"^] = 

for < A; < J — 1 and 

vL[^jZ-^]= v[A^jZ-^^+^'^]= v[{Cj^l^j+r + Dji;*^,)z-^^+'-^] = Cj{0)ej+r. 

So, Cj(0) / for j = 0, . . . , n iff {^j)'^^Q is a finite segment of MOP with 
respect to fL, which means that u = vL in P„. □ 
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Equality is true taking Dj = D*''^^ , no matter whether Dj has 
degree r — 1 or not. In what follows we will assume this convention for the 
polynomials Dj. 

Remark 2.5. The functional u has a finite segment of MOP of length (at 
least) one iff u[l] ^ 0. Therefore, Theorem 12.41 assures that the condition 
v[L] / is equivalent to the existence of a (unique) decomposition 

A = CoiJr + Dorr, CoGP., Do£Fr-l, (5) 

with Co(0) 7^ 0. However, Remark 12.31 says even more: no matter the 
value of v[L], there is always a unique decomposition like ([S]). The equality 
v[L] = Co{0)er implies that v[L] / is only responsible of Co(0) / 0. 

The above theorem has the following consequence for quasi-definite func- 
tional. 

Corollary 2.6. Letu, v be quasi- definite functionals with sequences of MOP 
{(fn), ii^n) respectively, and let L = P + = z~''' A with P a polynomial of 
degree r. Then, u = vL iff there exist polynomials C„ G Pj.; Dn G Pr-i with 
Cn{0) / such that 

Aipn = Cni^n+r + -DnC+r) « > 0, (6) 

or equivalently 

^fn = zD^i'n+r + C*V'*+r, n > 0. 

For convenience, in what follows we will use a matrix notation and we 
will adopt some definitions and conventions that will be used in the rest of 
the paper. If L is a hermitian Laurent polynomial of degree r, P and A are 
the polynomials given hy L = P + P^: = z~^A, P{0) E M. We denote by 
<f)j and ipj the j-th MOP with respect to the hermitian functionals u and v 
respectively. Also, 

Uj = (pj{0), bj = V'j(O), Cj = {ipj,ipj)u, Sj = {i;j,'ipj)y, 



^ zD* q J ' \D* c* 
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The matrices Sj and Tj, known as transfer matrices, permit us to write 
recurrence relations ([1]) and ([2]) for {<Pn) and (V'n) in the compact form 

^.='5,^,-i> ='?^-^j-i> (7) 

while the matrices Cj make possible to combine ([3]) and ^ into 

The structure of the matrices Cj is worth to be remarked. 
Definition 2.7. A polynomial matrix C = ^ C2 ) ' ^ satisfying 

C*'' = JCJ with J = ( ? n ) ^i^^ called a J- self- reciprocal matrix in P,.- 



1 0^ 

This is equivalent to state that C2 = C^'' and D2 = D^"^ . 

We denote by J,, the set of J-selfreciprocal matrices in P,. such that 
C2(0) 7^ and D2{0) = 0. These conditions mean that degCi = r and 
degl?! < r — 1, thus the general form of a polynomial matrix C G is 

C=f.n* n*)^ degC = r, degD<r-l, (8) 



^ C* ^ 

where here and below we assume that D* = D*^~^. 

Given a polynomial matrix C G Ir like ([8]) we will denote 



C 



C zD 

D* C* 



which is J-self-reciprocal too, but in general does not necessarily belong to 
because zD can have degree r. 

The determinant of a J-self-reciprocal matrix C in P^ is a self-reciprocal 
polynomial in ¥2r- When detC has degree 2r we will say that C is a regular 
J -self-reciprocal matrix. This is equivalent to detC(O) 7^ 0, which in case of 
C G Jr means simply C(0) / 0. We will denote by the subset of regular 
J-self-reciprocal matrices of J^- 

The next result about J-self-reciprocal matrices will be useful later on. 

Lemma 2.8. Let S = ( ^_ "'],T=(% ^) with a,beC. 

\za 1 J \zb 1 J 

(i) If\a\ 7^ 1, C G JJf.) the equation CT = SC defines a matrix C G Jr iff 

aC*{0) = bC{0)+D{0). 

In this case C e I?^ ^ |6| / 1, C G J'"^. 
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(ii) // 16| / 1, C G Jr, the equation CT = SC defines a matrix C £ iff 



aC(0) = bC*{0) - D*{0). 
In this case C £ fr""^ ^ \a\ ^ 1, C £ J^^. 
Proof. If \a\ 7^ 1 the equation CT = SC can be written as 

Let C £ Ir- Then A:" is a J-self-reciprocal matrix in P^, i.e., X = 

with X,Y £ Fr- Therefore, C is a polynomial matrix iff 1^(0) = 0, which 
yields the relation between a and b given in (i). In such a case Y = zY, 

Y £ Fr-i, and X*-{0) = C*{0) / 0, thus C = ^^t. ) £ Jr- Also, 

X(0) = C(0)(1 - - lap), hence C £ fr"""^^ ^ \b\ ^ 1, C £ J^'^. 

On the other hand, if |6| ^ 1 the equation CT = SC reads as 

1-|6|2 \a lj^\-b 1 

Suppose that C £ Ir- Then C is a J-self-reciprocal matrix in P^, hence 
C = ^y^r ^ with X,Y £ P,.. The relation between a and b given in (ii) 
is equivalent to y*'-(0) = 0, and also gives X*-{0) = C*(0) / 0, X{0) = 
C(0)(1 - |ap)/(l - |6|2) / 0, so C £ Ir and C £ fr""^ ^\a\^l,C £ fr""^. □ 

The goal of the rest of the section is to present a more economical and 
effective approach than the ones already existing in the literature (see for 
instance [HI [ISl [E]) to study the relation u = vL for any degree of L. 
This new point of view avoids the calculation of determinants and MOP 
related to u and v, requiring only the knowledge of the corresponding Schur 
parameters and the Laurent polynomial L. More precisely, we will charac- 
terize the relation u = vL through a matrix difference equation for the Schur 
parameters involving J-self-reciprocal matrices. 

The first step to formulate this new approach is to translate the relations 
between the MOP ((/?„) and (ipn) into relations between the corresponding 
Schur parameters. The following result will be useful for this purpose. 

Lemma 2.9. Let P, Q be relatively prime polynomials with degQ ^ degP. 
If the polynomial matrices 
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satisfy deg(M2 - A^2),deg(M4 - A^4) < degP, then 



Proof. MiP + M2Q = NiP + N2Q, thus (Mi - Ni)P = {N2 - Ma)^. Since 
gcd(P, Q) = 1, necessarily P divides M2 — N2, which impUes M2 — N2 = 
because deg(M2 — N2) < degP. Therefore Mi — Ni = too. Analogously 
M3 - N3 = M4 - Ni = 0. □ 

The next result is the matrix form of Theorem 12.41 together with a 
stronger result and some properties of the polynomial matrices Cj, including 
the first relations between the Schur parameters (0^) and (bn)- 

Theorem 2.10. Let u, v be quasi- definite in P„, Pn+r respectively and 
let L be a hermitian Laurent polynomial of degree r. Then, the following 
statements are equivalent: 

(i) u = vL in P^. 

(a) There exist Cq , . . . , C„ G II^^ such that 

A^j=Cj^j+r, j = 0,...,n. (9) 

(Hi) There exists Cn G Jr°^ such that 

A^n = Cn*n+r. (10) 

The matrices Cj are the only solutions of (0) in Sr, so Cq is determined by 

Co^'r = ^ I ) , CoeSr. (11) 

Besides, we have the relations 

CjTj+r = SjCj^i, j = l,...,n, (12) 

CjBj+r = AjCj^i, j = l,...,n, (13) 
detCj = Cj{0)A, j = 0, . . . , n. (14) 
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Proof. Bearing in mind Theorem 12.41 it is enough to prove (iii) =^ (ii) =^ 
(fni) . ([13]), (HI]). Suppose that only (iii) holds. Evaluating at z = we 
find a„A(0) = bn+rCn{0) + Dn{0) and C*(0) = A{0). Hence, Lemma[2S](i) 
assures the existence of C„_i G Sr°^ satisfying CnTn+r = ^nCn-i. Then, the 

equality ASn^n-l = A^n = Cn^n+r = CnTn+r^n+r~l = SnCn-l^n+r-l 

shows that A(^n~i = Cn-i'^n+r-i- Iterating this procedure we obtain (ii). 
Combining ^ and recurrence relations ([7]), 

A^j = Cj^j+r = CjTj+r^j+r-l, A^j = ASj^j-i = SjCj-i^ j+r-1- 

Therefore, CjT^+r^'j+r-i = SjCj^i'^j^r-i, or equivalently 



Taking into account that zipj, ip* are relatively prime and degCj = r, 
deg Dj < r — 1, relations (I12p and (jlSh follow from Lemma l2. 91 

To prove (fHl) notice that A = Cgipr + -Do^* = C^ip* + zD^ipr, hence we 
have the equality (Co — zDQ)ipr = (Cq — Dq)iP*.. Since ipr, ip^ relatively 
prime this implies Co{0)ipr = Cq — Do and Co(0)V'* = Co — zDq. So, 

Co(0)^ = Co(0)(Co^. + Dorr) = CoCo* - zDoD*^ = detCo. 



Besides, from ([12]) we find that detCj oc detCo for j = I, . . . ,n. Evaluating 
at z = we finally obtain detCj = (Cj(0)/Co(0)) det Co = Cj{0)A. □ 

The equivalence (i) 44> (iii) of the previous theorem means that the last 
condition {j = n) in ([3]) or ([!]) suffices for the equivalence in Theorem 12.41 

There exist also inverse relations between the finite segments of MOP 
ifj)^=o i'^j)'j=o- The polynomial matrix coefficients of these inverse 
relations are not independent of the polynomial matrix coefficients Cj of the 
direct relations. Indeed, both polynomial matrix coefficients are essentially 
adjoints of each other, understanding the adjoint of a 2 x 2 matrix M = 
as the matrix Adj(M) = ( "^f ). Thus, given a 2 x 2 



polynomial matrix M in P^, Adj(M) is a 2 x 2 polynomial matrix in P,. 
satisfying 

Adj(M)M = (detM)/, 
where / is the identity matrix of the same size as M. 

Theorem 2.11. If u, v are quasi- definite m P„, Pn+r respectively, the fol- 
lowing statements are equivalent: 
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(i) u = vL in P„ for some hermitian Laurent polynomial L of degree r. 
(a) There exist Xr, . . . , Xn S ^^r^ such that 

^j+r = Xj^j, j=r,...,n. (15) 

(Hi) There exists Xn £ such that 

"^n+r = Xn^n- (16) 

The matrices Xj are the only solutions of (f73|) in Jr, so Xr is determined by 

Xr^r = '^2r, G Jr- (17) 

Besides, we have the relations 

Tj+rXj^i = XjSj, j = r + 1,. . . ,n. (18) 

Bj+r'^j-i = ^j-^jj j = r + 1, . . . ,n. (19) 
detXj(xA, j = r,...,n. (20) 

CjX,=A(^^^ j = r,...,n. (21) 

X, = -^Adi{q), j = r,---,n. (22) 

Proof. If u = vL in P„, Theorem 12.101 assures the existence of Cj G Jr*^^ 
such that A^j = Cj^j-^r for j = 0, • • • ,n. Multiplying this identity on the 
left by Adj(Cj) and taking into account (fT^ we find that ^'j+r = '^j^j 
for j = 0,...,n, where Xj = Adj(Cj)/Cj(0) G J^^. Then, (fTSj) . (fT9l). (pO]) 
and (j2ip are a direct consequence of (jl2p . (jl3p and (|14p . The uniqueness 
of Xj £ Sr for j > r follows from Corollary 12.21 and the fact that (jl5p is 
equivalent to V'i+r = ^j^j + ^Vj) where Xj G Pj-, S Pr-i are the 

polynomials appearing in Xj = ^ ^ . 

It only remains to prove (iii) (i). Multiplying ()16p on the left by 
C„ = Adj(^„) G Sr°^ we obtain A^n = Cn^n+r where A = detA'„ is a 
self-reciprocal polynomial of degree 2r. This proves that u = vAz~^ due to 
Theorem EH □ 
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Concerning the polynomial matrix coefficients Xj G Sl^^ of the inverse 
relations, when it is necessary we will use the explicit notation 

= ( ^* X* ) ' ^"^^ " ^"^^ = ^ - 1' 

so that ^j+r = ^j^j is equivalent to ipj+r = ^jfj + ^j^j- This shows 
that Xj is monic. Besides, from (j22p we have the relations Xj = C*/Cj{0), 
Yj = -Dj/C,{0). 

The proof of the previous theorem shows that, when u = vL in for 
some hermitian Laurent polynomial L of degree r, 

"Sj+r = Xj^j, G Jr, j = 0, . . . , n, (23) 

and not only for j > r. Indeed, the proof of the theorem implies that (I23p 
has solutions Xj G Sr^^ for j < r too. The only difference is that, contrary 
to j > r, (f23j) does not determine univocally for j < r, as Corollary 
12.21 points out. The reason is that «B = {z^^jl^^o U {z'^V'jlllo is hnearly 
independent for j > r, but not for j < r. Actually, when j < r, Lemma |2. II 
shows that rank(!B) = j + r + 1, so the solutions Xj of (I23p form an affine 
subspace of dimension r — j. 

Among the solutions of ()23p for j < r there is a choice of special inter- 
est: similar arguments to those at the beginning of the proof of Theorem 
12.101 show that Lemma 12.81 (i), together with (I17p . assures that (|18p can be 
extended in a unique way to j = 1, . . . , r, giving rise to particular solutions 
Xq, . . . ,Xr-i G Jr*^^ of (f23]) . The choice of Xj determined by the extension 
of (fT8|) has the particularity that det Xj is independent of j up to numerical 
factors. Indeed, this property characterizes such a particular choice be- 
cause different solutions of (j23p can not have proportional determinants: let 
A'W, G Sr be such that = X^f'^^j. Then, (det A'W)^>j = C^'^^'^j+r 
with C^'^) = AdiiX^''^). If det = Adet A'(^), A G R*, Lemma ED assures 
that = XCW, thus detA:'^^) = X'^detX^^ which implies A = 1, so 

x^^^=xw. 

Properties (|12p and (jlSp are the cornerstone of the main objective of 
this section: a new characterization of the relation u = vL in terms of a 
recurrence for the corresponding Schur parameters. Like in the previous 
characterizations, the J-self-reciprocal matrices play an important role, but 
now only one MOP of u and v enters in the equivalence, and it appears only 
in the initial condition for the recurrence. The direct and inverse relations 
between the MOP of u and v lead to different characterizations, depending 
on whether the hermitian Laurent polynomial L is fixed or not. Indeed, L 
appears explicitly only in the initial condition for the direct characterization. 
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Theorem 2.12. Let u, v be quasi- definite in ¥n, Fn+r respectively and 
consider an index m G {0, . . . , n}. 

(i) Given a hermitian Laurent polynomial L of degree r, u = vL in ¥n iff 
there exist Cm G and Cm+i, • • • ,C„ G such that 

C-rri^m+r = ^^m, (Direct Initial Condition) 

CjTj^r = SjCj^i, j = m + I, . . . ,n. (Direct Recurrence) 

Moreover, A^j = Cj'^j^r, Cj G Sl^^ and detCj oc A for j = m, . . . ,n. 

(a) There is a hermitian Laurent polynomial L of degree r such that u = vL 
in ¥n iff there exist Xm G and Xm+i, ■ ■ ■ , Xn G Jr such that 

'^m^m = ^m+rj (Inverse Initial Condition) 

Tj+r'^j-i = ^j'Sj, j = m + 1, . . . ,n. (Inverse Recurrence) 

Moreover, ^j+r = ^j^jj G Sl^^ and det Xj oc A for j = m, . . . ,n. 

Proof. We will prove only (i), the proof of (ii) being similar. In view of 
Theorem 12.101 it suffices to show that Direct Initial Condition and Di- 
rect Recurrence imply A^j = Cj'^j+r, j = m,...,n and Cj G Jr°^, j = 
m + 1, . . . ,n when Cm G Sl^^. Direct Recurrence yields (1 — |6j+,.p) detCj = 
(1 — |ajp) detCj_i, thus Cm G Sr°^ implies Cj G Jr°^ for j = m + 1, . . . , n. 
Also, Direct Recurrence and Direct Initial Condition combined with recur- 
rence relations ([7]) lead to A^j = ASj • • • Sm+i^m = "Sj • • • Sm+iCm'^m+r = 

CjTj+r ■ ■ ■ %n+r+l^m+r = Cj'^j+r foi j = m, . . . ,n. □ 

Some special cases of the above theorem will be of interest for us. We 
will summarize them. 

Theorem 2.13. Let u, v be quasi-definite in P„, IPn+r respectively. 

Direct characterization Given a hermitian Laurent polynomial L of degree 
r, the following statements are equivalent: 

(i) u = vL in P„. 

(ii) There exist Cq G and Ci, . . . ,Cn & such that 




(Initial Condition D) 



CjTj+r = SjCj-i, j = 1, . . . ,n. (Recurrence D) 
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Inverse characterization The following statements are equivalent: 

(i) u = vL in P„ for some hermitian Laurent polynomial L of degree r. 

(a) There exist G Jr'^^ and Xy-\-i, . . . , Xn G Jr such that 

Xr^r = *2r, (Initial Condition II) 

Tj^r^j~i = ^j'Sj, j = r + 1, . . . , n. (Recurrence II) 
( Hi ) There exist Xq € Ir^^ and Xi, . . . , Xn G Jr such that 

;fo ^ I ^ =^r, (Initial Condition 12) 

Tj^r-^j-i = '^j'Sj, J = 1; • • ■ 5 (Recurrence 12) 

The difference between the inverse characterizations II and 12 is that the 
initial condition determines univocally the initial matrix Xr for II but not 
the initial matrix Xq for 12, thus there is a freedom in such initial matrix 
for 12. We will go back to this point later on. 

Theorems 12.101 12.111 12.121 and 12.131 have an obvious generalization to the 
quasi-definite case. 

Theorem 12.131 shows that the regularity of Cj, j ^ 0, and Xj, j ^ r, is 
a superfluous condition in statement (ii) of Theorems 12.101 and 12.111 respec- 
tively. Remember that the regularity of Cq is equivalent to v[L] ^ 0. On 
the other hand, the regularity conditions for Xj in Theorems 12.111 and 12.131 
can be completely avoided if we do not fix the degree of L. In other words, 
if Xj G Jr \ Jr'^^ then u = vL in P„ too, but degL < r, as follows from the 
following proposition. 

Proposition 2.14. If^j+r = ^j^j with Xj G Sr-Xfr""^, then^j^r-i = ^j^j 
with Xj G Jr-i- 

Proof Suppose ^j+r = ^j^j, G Jr with Xj(0) = 0. Then hj+r = ^i(O) 
and Xj = zXj with Xj monic of degree r — 1. Thus we can write 

From ^j+r = Tj+T^ i+T-i we get ^'j+r-i = Xj^j where 

_ 1 / X, - bj+rY; Z-\Yj - bj+rX*) ■ 

1 - \bj+r\^ V Z{Y* - bj+rXj) X* - bj+rYj 

Since Yj{0) - bj+rX*{0) = and X*{0) - bj+rYj (0) = 1 - |6j+rP / we 
conclude that Xj G Jr-i. □ 



X'i 
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3 Direct and inverse problems 



In the previous section, given two hermitian linear functionals n, v and 
a hermitian Laurent polynomial L, we have studied the relation u = vL 
obtaining characterizations in terms of linear relations with polynomial co- 
efficients between the corresponding MOP, as well as in terms of a matrix 
difference equation between the related Schur parameters. In this section 
we will use these results to answer the following question: Which conditions 
ensure the quasi-definiteness of n = vL or v once we know that the other 
functional is quasi-definite? 

Indeed we will answer this question in the more general context of quasi- 
definite functionals in some subspace P„: we will try to know the minimum 
length of the finite segments of MOP for one of the functionals assuming 
that the other functional has a finite segment of MOP with a given length. 
Like in the previous section, the main goal is to develop techniques for this 
problem based almost exclusively on the knowledge of the Schur parameters. 

The new results will seem quite similar to those of the previous section, 
however they provide new information: in the previous section we assumed 
that u and v had finite segments of MOP of certain length and we asked 
about a characterization of the relation u = vL in some subspace P^; now 
we will consider the relation u = vL as a data and we will ask about the 
length of the finite segments of MOP. 

3.1 Direct problem 

The direct problem refers to the case where we suppose that a hermitian 
functional v with a finite segment of MOP {ipj)JLQ and a hermitian polyno- 
mial L of degree r are given. Then, we will try to obtain information about 
the functional u = vL and its finite segments of MOP (v7j)"=o- ^^^^ 
result is essentially a reinterpretation of relation ([9]) . 

Theorem 3.1. Let v be quasi-definite ifi IPy^_^^ ciTid let L he CL hevuiitidTi 
Laurent polynomial of degree r. Then, u = vL is quasi- definite in P„ iff 
there exists Cj £ ll^^ such that A divides Cj"^j+r for j = 0, . . . ,n. 

Besides, det Cj oc A and there is a unique choice of Cj such that Cj (0) = 
^(0). For such a choice the finite segment of MOP with respect to u is given 
by A<^>j = Cj^j+r for j = {),..., n. 

Proof. First of all notice that, no matter the value of \j G C*, ^ divides 
Cj^'j+r- iff it divides Cj^'j+r with Cj = ^ 'q' Cj, and Cj G Jr'^^ iff Cj G Jr'^^. 
Therefore, we can suppose without loss of generality that C*(0) = ^4(0). 
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Then, the divisibihty condition is equivalent to A^j = Cj^j^r with (pj a 
monic polynomial of degree j which, for the moment, has no relation with 
u. Taking into account Theorem 12.101 to prove the result we only need 
to see that (pj is the j-th MOP with respect to u. This follows from the 
orthogonality conditions of 'i/'j+r with respect to v, which give 

u[p,z->']= v[{C,i,,+r + Djrj+r>-^''^"^] =0, < A; < i - 1, 

The rest of the theorem is a consequence of Theorem 12.101 □ 

The above results allow us to obtain a necessary and sufficient condition 
for the quasi-definiteness of the functional u = vL in terms of determinants 
involving the MOP of v. 

Proposition 3.2. Let v be quasi-definite ifi IP^_|_^. and let L be a hevmitian 
Laurent polynomial of degree r. Then, u = vL is quasi-definite in P„ iff 
det M("^) / /or m = 0, . . . , n + 1, where M^'") = (M^"'^)^^.^^ is the square 
matrix of order 2r given by 

^ 1 {z^-'^Pn.+rt{a j = h...r, , = 1, . . . , 2r, 

\{z^~'^'rm+rf'{a j = r + l,...,2r, 

with <^i, . . . , ^2r the roots of A counting the multiplicity and U the number of 
roots Qj, j <i, such that C,j = Q. 

Proof. By Theorem 13.11 to decide the quasi-definiteness of u = vL in P.„, we 
simply have to analyze the existence of unique polynomials C^, -D^ with 
degCm = r, degDm < 1, (7^(0) 7^ 0, C^{^) = A(0), such that A divides 

Cmlpm+r + -DmV'm+r for ?n = 0, . . . , n. 

Let us write Cm{z) = Ylk=o^rn,kz'' and Dmiz) = Ylk=odm,kz'' ■ The 
condition C^(0) = A{0) only means that Cm,r is the leading coefficient 
of A. Then, the existence of unique polynomials Cm, Dm is equivalent 
to the existence and uniqueness of the 2r coefficients Cm,o, ■ ■ ■ ,Cm,,r-i and 
dmfl, • • • , dm,r-i, while the condition Cm(0) 7^ becomes Cm,o / 0. 

If Ci, • • • , C2r denote the 2r roots of the polynomial A counting the mul- 
tiplicity and li is the number of roots Q such that ("j = Q for j < i, the 
divisibility condition is equivalent to the system 

{Cmtpm+r)^''{Q) + (Ci) =0, i = 1, . . . 2r. 
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This system has a unique solution in Cm.k, dm,k, k = 0, . . . ,r — 1, exactly 
when det M^™) / 0. 

It remains to translate the condition Cmfl 7^ 0. The solution for Cmfl 
is proportional to the determinant of a matrix M obtained substituting in 
M("^) the first column (^S'+^lCi))?!! by (z^^2:+^(Ci))?Ii. Since © and © 
imply that span{z-''+Vm+r, ^•'V'm+r} = span{z^ ipm+r+i, , we see 

that det M vanishes at the same time than det M^™'^^) . Hence, Cm,o 7^ is 
equivalent to det M(™+i) ^ 0. □ 

The condition given by the above proposition is theoretically interesting 
but in practice it is not manageable, specially for polynomial perturbations 
of high degree r due to the need to evaluate determinants of 2r x 2r matrices. 
Even in case of low degree r, the practical application of the previous result 
needs the construction of the MOP ipj and the evaluation at some points of 
these MOP and their derivatives. 

When r = 1 the self-reciprocal polynomial A has two roots Cij C2 such 
that C2 = 1/Ci or Ci, (^2 G T, (i (2- Obviously, when v is positive def- 
inite and = 1/Ci the functional vL is positive definite too. However, 
in general, v quasi-definite in Wji-^j- implies vL quo-si-dcfinite in IP^^ iff (sgg 
pa a El [1]) Km{Ci, l/C2) / for m = l,...,n + 1, where Km{z,w) = 
^Y=o^J^'4'j{z)ipj{'w) is the m-th kernel associated with the MOP (ipj). 

Nevertheless, it is naive to think that the general situation can be solved 
by factoring the polynomial A. Consider for instance a positive definite 
functional v and let A{z) oc {z — Ci)iz — C2) with C11C2 G T, (^1 / ^2) 
satisfying Km{Ci, I/C2) = for some m. Then vL is not quasi-definite but 
vL'^ is positive definite. 

A more practical characterization of the quasi-definiteness of n = vL, 
which avoids the construction of the MOP of v and does not need the cal- 
culation of determinants, is given in terms of the recurrence for the Schur 
parameters. 

Theorem 3.3. Let v be quasi-definite ifi IP^^y-. and let L he a hevmitian 
Laurent polynomial of degree r. Then, u = vL is quasi- definite in P„ iff 
there exist oi, . . . , a„ G C and Cq, . . . ,Cn G Jr°^ such that 

Co^r = A(^\y (24) 

CjTj+r = SjCj-i, j = 1,... ,n. (25) 

Besides, A^j = Cj'^j+r, detCj oc ^4, j = 0, . . . , n, and aj = (pj{0) G C \ T, 
j = l,...,n. 
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Proof. In view of Theorem 12.131 we only need to prove that u is quasi- 
definite in P„ when dMD and 1^ hold. Define = Sj---Si(\). Then, 



23])) (|25p and the recurrence relation for {^j)"'^Q yield for j = 0, . . . ,n, 



A<^j = ASj ■ ■ ■ Si [ ) = Sj ■ ■ ■ SiCo^r = CjTj+r ■ ■ ■ Tr+l'^r = Cj^j+r- 



Therefore, Theorem 13.11 shows that u is quasi-definite in P„. □ 

The above results yield a direct relation between the Schur parameters 
oi u = vL and v, which can be obtained setting z = in the equivalent 
version CjBj+r = AjCj^i of ([25]) and using C*(0) = ^(0). 

Corollary 3.4. If a is the leading coefficient of A, thej-th Schur parameter 
Qj of u = vL can be obtained from the j + r-th Schur parameter bj+r of v by 



_ abj^r-D*.m 



(26) 



Theorem 13.31 and Corollary 13.41 provide an algorithm to obtain the Schur 
parameters (oj) of u = vL from the Schur parameters {bj) of v. 

Algorithm D 

• Determination of Co G Jr from initial condition ()24p and A. 

• For j = 1,2,... 

• While Cj_i(0) 7^ 0, calculation of Cj from (j26p and ^j-i- 

• Determination of Cj G from recurrence ([25]) and Oj, 6j+r, Cj-i. 

The fact that the j-th step of the above algorithm actually gives a matrix 
Cj G Jr is a consequence of Lemma 12.81 (ii) and the equivalence between 
Cj_i(0) / and Cj-i G fr^ when Cj-i G J,-. 

In short, the fact that Algorithm D works from j = 1 to j = n will be 
called the n-consistence of recurrence (|25p . We will say that the recurrence 
is consistent if it works for any j > 1. Of course, this is an abuse of language 
because the consistence depends, not only on recurrence ([25]) . but also on 
initial condition ()24p . 

The consistence relies on the fact that Cj(0) 7^ at each step. Suppose 
that the recurrence fails at the (n+ l)-th step, i.e., it is ?i-consistent and not 
(n + l)-consistent. Then C„_i(0) 7^ and C„(0) = 0, that is, C„_i G J: 



reg 
r 
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but Cn G Jr\Jr°^- Recurrence ()25p shows that this is equivalent to |an-i| 7^ 1 
and |a„| = 1. So, the n-consistence condition can be written as \aj\ ^ 1 for 
j = 1, . . . ,n — 1, which means that u = vL has a finite segment of MOP of 
length n, i.e., it is quasi-definite in Pn-i. 

Contrary to Theorem 13.21 Algorithm D only requires the knowledge of 
the Schur parameters of v and a single MOP tpr with the same degree r 
as the polynomial perturbation L. Furthermore, this algorithm makes the 
calculation of determinants completely unnecessary. As an example, we will 
develop explicitly Algorithm D for r = 1. 



3.1.1 The case r = 1 

Consider a hermitian functional v with MOP (ipj) and a hermitian Laurent 
polynomial L of degree 1. We can write L = P + P^^, P{z) = az + (3, a & C* , 
/? G M, so A{z) = zL{z) = az^ + 2(3z + a. The MOP {tpj) of the modified 
functional u = vL, if they exist, are given by 

Aipj = {az + Cj)^n+i + djip*+i, 
for some Cj S M, dj £ C. This relation and its reversed can be combined in 



Also, recurrence (j25]) becomes 



az + Cj dj 
djZ a + CjZ 



Cj„i + dj-ittj = Cj + djbj+i, 
aaj = Cjbj + dj, 
Cj-iaj + dj-i = a6j+i. 



which can be written as 



abj+i - dj^i / 1 bj+i \ f Cj \ _ f Cj-i + dj-iQj 



On the other hand, initial condition (|24p is equivalent to az"^ + 2(3z + a = 
{z + ci)(z + 6i) + doihz + 1), i.e., 

This provides unique cq, do for any P and any possible value of 6i S C \ T. 
Finally, Algorithm D can be explicitly formulated in the following way: 
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• Calculation of cq, do from P, bi using (p8]) . 



• For j = 1, 2, ... , while Cj_i 7^ 0, calculation of aj, Cj, dj from 6j+i, 
Cj„i, using dSZI). 

This algorithm provides the Schur parameters of u = vL and informs us 
about its quasi-definiteness: the maximum subspace P.„ where u is quasi- 
definite is given by the first index n of inconsistency of the algorithm. 

We can think in reducing the general problem to the case r = 1 by factor- 
ing the polynomial A. Suppose that A = A1A2, degAi = 2ri, deg = 2r2, 
with Ai self-reciprocal, and denote by Cj^^ , cj^^ the J-self-reciprocal matrices 
associated with the direct problem w = vAiz~^"^, u = wA2Z~^'^ respectively. 
If Uj are the transfer matrices for the functional w with MOP (^j) and 

Ej = ^1*^1 then AiEj = Cj^^^j^n and A2^j = cj^^Hj+r2- This implies 

the equality A^j = cf'^C^^2^ ^i+rj so Cj = cf'^C^^^. . However, this does 
not always reduce a direct problem to simpler ones because the length of 
the finite segments of MOP for w can be not big enough to get the actual 
relations between all the MOP of u and v. 



3.2 Inverse problem 

In this subsection we will study a problem which can be consider as the in- 
verse of that one of the previous section. More precisely, given an hermitian 
functional u with a finite segment of MOP {ipj)^^Q and a hermitian Lau- 
rent polynomial L of degree r, we will try to obtain information about the 
hermitian solutions v u = vL and their finite segments of MOP (V'j)jLo- 
First of all we will clarify the structure of the set {v hermitian : u = vL} . 
The equation u = vL is equivalent to ^[2;"] = v[z'^L], n > 0, which, denoting 
Hn = u[z'^], rUn = f [z"] and L{z) = Yl^j=-r otj^-' O-j = aj, becomes 

r 

= ^ ajuin+j, n > 0. (29) 

j=-r 

The first equation (n = 0) 

r 

//Q = 2 Re ^2 (30) 

j=0 

is simply a constraint between the first r + 1 moments uiq, . . . ,mr of v. 
The rest of the equations determine the moments m„, n > r. Since any 
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hermitian solution v is determined by its moments m„, > 0, the general 
solution depends on 2r real independent parameters obtained establishing 
in the set {mo, mi, . . . , m^}, mg G M, mi, . . . , m^ G C, the constraint (j30p . 

There is another way to describe the set of hermitian solutions v starting 
from a particular one vq. Then the hermitian solutions are those functionals 
with the form u = + A, where A is any hermitian functional satisfying 
AL = 0, i.e., 

^ = E E - < ^ c, ii4f = M« if = i/c., 

Cj, i = being the roots of ^ = z~'^L and the multiplicity of 

C,i. Again we see that the hermitian solutions are parametrized by 2r real 
parameters: the independent real and imaginary parts of the coefficients 
. Furthermore, this approach shows that the inverse problem is related 
to the study of the influence of Dirac's deltas and their derivatives on the 
quasi-definiteness and the MOP of a hermitian functional. 

For convenience we will denote by H^. [u) the set of hermitian functionals 
V which are solutions oi u = Lv for some hermitian Laurent polynomial L 
of degree r. The main result of this section characterizes the functionals of 
Hr{u) which are quasi-deflnite in some subspace Pm- 

Theorem 3.5. Let u be quasi- definite in P„. 

(i) If n > r, there is a (unique up to factors) solution v € Hr{u) quasi- 
definite in Fn+r for each 6i, . . . , 62r G C \ T, 62r+l, • • • , ^n+r G C, 
Xj., . . . , Xn G Jr^^ such that 

Xr^r = ^2r, *2r = " " " ^i J ^ , (31) 

Tjj^r^j-i = j = r + 1, . . . , n. (32) 

The relation between v and bj, Xj is that ^j+r = ^j^j provides the 
j + r-th MOP of V for j = r, . . . ,n, and bj € C \ T, j = 1 , . . . , n + r , are 
the first n+r Schur parameters of v. Besides, u = vL with detXj oc A, 
j = r,...,n. 

(ii) There is a (unique up to factors) solution v S Hr{u) quasi-definite in 

¥n+r for each bi, . . . ,br £ C\T , br+l, ... , bn+r £ C, Xq, . . . ,Xn £ Jr°^ 

such that 

^o(})=^'r, ^r=Tr---Tj\Y (33) 
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Tj+rXj-i = XjSj, j = 1, . . . , n. (34) 

The relation between v and bj, Xj is that ^j+r = ^j^j provides the 
j + r-th MOP ofv for j = 0, . . . , n, and bj € C\T, j = 1, . . . , n + r, are 
the first n+r Schur parameters ofv. Besides, u = vL with detXj oc A, 
j = 0, . . . ,n. 

Proof. We will prove only (i), the proof of (ii) is similar. Bearing in mind 
Theorems 12.111 and 12.131 we only need to show that (j3ip and (j32p imply that 
Xj^j gives for j = r, . . . ,n the j + r-th MOP of a unique v G Hr{u) whose 
first n + r Schur parameters are bj, j = 1, . . . ,n + r. 

Let us define = 7^- • • • 71 ^ | ^ Since Xj £ Jr°^, j = r, . . . ,n, recur- 
rence (j32p implies that \bj\ ^ 1, not only for j = l,...,2r, but also for 
j = 2r + 1, . . . , n + r. Therefore, (''/'i)j=o ^ finite segment of MOP with 
respect to some hermitian functional v. 

From ([3T]) . ([32]) and the recurrence relation for ($j)"=Q we obtain for 
j = r,...,n, 

Hence, Theorem 12.111 proves that u = vL in IPy^ for some hermitian Laurent 
polynomial L of degree r. Multiplying L by a real factor we can get a 
hermitian Laurent polynomial L of degree r such that u = vL in P„. 

The equality n = vL in P„, as well as the fact that {ipj)^^Q is a finite 
segment of MOP for v, only depends on the first n + r + 1 moments v\z^], 
J = 0, . . . , n + r, of {). Let us define a new hermitian functional u fixing its 
moments rrij = v[z^] by rrij = v\z^] for j < n + r, and ruj given by (I29p for 
J > n + r + 1 . Then u is a solution of n = vL, has (^j )"^q as a finite segment 
of MOP and its first n + r Schur parameters are ipj{0) = bj, j = 1, . . . ,n + r. 

Finally, the first n + r Schur parameters of a functional v determine its 
finite segment of MOP of length n+r+1 and, thus, its first n + r + 1 moments 
up to a common factor. Requiring also u = vL for a given hermitian Laurent 
polynomial of degree r fixes the rest of the moments up to the common 
factor due to ()29p . Therefore, the conditions of (i) define a unique hermitian 
functional v up to factors because L is determined up to real factors by 
detXj. □ 

We have the following relation between the Schur parameters of u and 
V £ Hr{u). To prove it simply choose z = in the equivalent version 
Bj+rXj-i = XjAj of ([M]) and use that Xj is monic, i.e., ^*(0) = 1. 
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Corollary 3.6. The j + r-th Schur parameter bj^r of v £ Hr{u) can be 
obtained from the j-th Schur parameter Uj of u by 



a,- - Y* . (0) 
bj+r = ' • (35) 

X,_i(0) 

Theorem 13.51 and Corollary 13.61 provide algorithms generating the solu- 
tions of the inverse problem which are quasi-definite in some subspace Pm- 
The algorithms are based on the consistence of recurrence (|32|) or (j34|) , what 
can be defined in a similar way to the case of Algorithm D. We have several 
possibilities depending of the initial data. 

If we know that L has degree r but not its explicit form, we can proceed 
in the following ways, depending whether we are interested in the solutions 
which are quasi-definite (at least) in or P^. 

Algorithm II 

• Choice of ^'2r, i-e., of 6i, . . . , 62r G C \ T. 

• Determination of G Jr from initial condition pip and <!>,., ^2r- 

• For j = r + 1, r + 2, . . . 

• While Xj-i{0) 7^ 0, calculation of bj^r from ()35|) and Cj, -^j-i. 

• Determination of Xj G from recurrence (|32p and Oj, 6j+r, ^j-i- 

Algorithm 12 

• Choice of ^'r, i.e., of 6i, . . . , 6,. G C \ T. 

• Choice of a solution Xq G Sl^^ of initial condition ([33]) using i-e., 
choice of a monic polynomial Xq of degree r with Xq{0) / and 
determination of Yq = ipr — Xq . 

• For j = 1,2,... 

• While Xj_i(0) 7^ 0, calculation of bj^r from ([35|) and Oj, ^j-i- 

• Determination of G Jr from recurrence (|34p and Oj, 6j+r, '^j-i- 

For any of these two algorithms we recover the polynomial perturbation 
through A oc det Xj . 

On the contrary, if we know explicitly the hermitian polynomial L of 
degree r, we have the following scheme to find the solutions v oi u = vL 
which are quasi-definite in P^. 
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Algorithm 13 



• Choice of i.e., of 61, . . . , &r G C \ T. 

• Determination of = '^^ from initial condition (j24p and ^• 

. For J = 1,2,..., 

• While Xj_i{0) 7^ 0, calculation of bj^r from (j35|) and a^, ^^j-i. 

• Determination of A'j G from recurrence p4p and aj, bj^r, ^j-i- 

We can assure that any step of the above algorithms generates a matrix 
Xj £ due to Lemma |2. 81 (ii) and the fact that Xj_i(0) 7^ is equivalent 
to Xj-i G Jr*^^ when Xj^i G Sr- 

The n-consistence of the above algorithms, which means that they work 
for j < n, is equivalent to the existence of a finite segment of MOP of length 
n + r for the corresponding solution of n = vL. Such n-consistence can be 
written as Xj{0) 7^ 0, j < n — 1, which holds iff \bj\ ^1, j<n + r — 1. 

Comparing the above algorithms we see that the arbitrariness in the 
parameters br+i, . . . , is equivalent to the arbitrariness of the polynomial 
modification L of degree r. This means that any of the infinitely many 
solutions Xq G Jr^^ of Xq^q = should be determined by det ^Yq, a result 
which is proved in the next proposition. 



Proposition 3.7. Given bi, ... ,br £ C\T and a self-reciprocal polynomial 



A of degree 2r, there exist a unique solution Xq G Sr°^ of Xq i | ) = "^r, 



= 7^. • • • 7i ^ I J , such that det Xq oc A. 

Proof. Given each solution of ^%'o*^*o = with the form 

^0 = (^^, ) , degXo = r, degFo < r - 1, 

is determined by a monic polynomial Xq because Yq = tpr — Xq. Then 
Xq G Jr'^ iff Xq{0) / 0. Therefore, det A'o = X^Vr + ^o^^ - Ai^r- Hence, 
if A is a self-reciprocal polynomial of degree 2r and A G M, 

det Xq = XA ^ XA + i)r^; = X^^r + ^oV'r ■ (36) 

From Remark 12.31 we know that \A + ipr'^* ~ Cipr + Di^r fo'^ some 
polynomials C G Pr, D G Pr-i- Since + V'rV'r is self-reciprocal in P2r, 
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(C - zD*)i;r = {C* - D)ip;, soC*-D = c^r and C - zD* = op; for some 
c G M. Then, the identity 

XA + iPrrr = iC- ^rr)A + {D + |Vr)€ = ^(^ + zD*)^r + ^{0* + D^r 

proves that (f36l) holds with Xq = \{C* + D), thus Xq satisfies det ^Yq = XA 
with such a choice. Furthermore, Xq is monic of degree r iff Xq{0) = 1, 
which (j36]) shows that corresponds to A = Xq{0)/A{0). 

Now, let A'o,*^b G Jr°^ be such that A'o^o = ^o^o = Assume that 
det Xq = X det A'o for some A G M. Using an obvious notation, this means 
that XQipr + (-'^o — V'r)V'* = M-^oi^r + (-'^o " V'r)V'r )• The uniqueness of the 
polynomials C, D in Remark 12.31 ensures that Xq = XXq and Xq — ipr = 
X{Xq — ipr), which implies that Xq = Xq. □ 

The previous results show that the solutions of the inverse problem are 
parametrized by their first r or 2r Schur parameters, depending on whether 
we fix the polynomial perturbation or only its degree. Of course, such a 
parametrization works only for the solutions which are quasi-definite (at 
least) in and P2r respectively. Each of these solutions will have a finite 
segment of MOP of maximum length determined by the consistence level of 
the corresponding algorithm. 



3.2.1 The case r = 1 

As an example of the previous discussion we will analize the particular case 
of the inverse problem corresponding to a hermitian Laurent polynomial 
perturbation L of degree 1. So, we consider the MOP {(pj) with respect to 
a hermitian linear functional u and we define the monic polynomials (ipj) 

ipj+i = {z + Xj)ipj + yjipj, j > 0, (37) 

with iPq{z) = 1 and xj, yj G C. The polynomials {ipj) are the only candidates 
to be MOP of a solution v oi u = vL. 
We can write ()37p in a matrix form as 

^,+1 = ^,$,, ^, = f" + !^' ,1'.), 3>0, (38) 



and (1331) becomes 



yjz 1 + XjZ 



Xj-i + Ijj-ibj+i = Xj + yjUj, 

bj+i = XjUj + yj, (39) 
Xj-ibj+i + Uj-i = aj, 
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or equivalently, 



\aj ^ J \ yjj V ^+1 J 

So, Algorithm 12 reads as follows: 

• Choice of 6i € C \ T and xq G C* which determines yo = bi — xq. 

• For j = 1, 2, ... , while Xj-i ^ 0, calculation of 6j+i, Xj, yj from aj, 
Xj-i, yj-i using (001). 

For any choice of xq we can recover the polynomial perturbation through 
A oc det = xqz^ + (1 + |xoP — lyoP)^; + Xq. According to Proposition 13.71 
given bi, each choice of xq in the previous algorithm provides a solution of 
the inverse problem corresponding to a different polynomial perturbation. 
These solutions have well defined MOP Voi ^i) so the algorithm provides 
all the solutions of the inverse problem which are quasi-definite at least in 
Pi. The maximum length of the finite segments of MOP for a particular 
solution is equal to the consistence level of the algorithm starting with the 
values bi and xq defining such solution. 

It is remarkable that, when r = 1, the consistence of Algorithm 12 is 
equivalent to the compatibility of ([Ml) , i.e., any solution of (f39]) for j < n 
starting with xq ^ necessarily satisfies Xj ^ for j < n — 1. We can see 
this by induction: if (I39p has a solution for j < n + 1, then Xn-i 7^ due to 
the induction hypothesis, so x„ = would give bn+i G T according to ([M]) : 
on the other hand, setting = in ([39l) for j = n, n + 1 we get y„ = bn+i 
and yn = 0^+1, which is a contradiction because 0^+1 ^ T. 

3.2.2 An example of the inverse problem 

As an application, we will solve the inverse problem for an arbitrary her- 
mitian polynomial perturbation L of degree 1, when u is the functional 
associated with the Lebesgue measure on the unit circle 

More precisely, we will characterize the quasi-definite solutions v G Hi[u). 
Indeed, we will do something more than this because our methods permits 
us to characterize all the solutions v G Hi[u) which are quasi-definite at 
least in Pi, providing also the maximum subspace Pm where each of such 
solutions is quasi-definite. 
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Bearing in mind the comments at the beginning of Section [3^21 and taking 
into account the possibiUties for the roots of a self-reciprocal polynomial A 
of degree 2, this is equivalent to the analysis of functionals v with the form 



(a) vo + M6{z - C) + M5{z - 1/C), C e 10, M £ C, 

(b) vo + Mi6{z - C) + M26'{z - C), C G T, Mi E M, 

(c) t;o + Mid(z-Ci) + M25(z-C2), C1/C2, CiGT, M, G M, 



where vq is a particular solution of the inverse problem. In case (a) we 
can take vq as a multiple of the functional associated with the measure 
dm{z)/\z — and then (a) is known as the Geronimus transformation 
of the Lebesgue measure. The Geronimus transformation of an arbitrary 
positive measure on the unit circle has been studied in [HI [9], while a more 
general laurent polynomial transformation has been analyzed in [241 [4j \W\ 
[TT|[5]. Our approach permits us to deal with the above three transformations 
simultaneously. 

The functional u is positive definite with MOP V3n.(-z) = -2", n > 0, and 
Schur parameters a„ = 0, n > 1, so that ([39]) becomes 



Following Algorithm 12, every choice of 61 E C \ T and xq E C* determines 
Vo = bi — Xq providing initial conditions for the above recurrence. Each of 
such initial conditions is associated with a different solution of the inverse 
problem we are considering, and this solution is quasi-definite exactly when 
the related initial conditions make ()4ip compatible for every n E N, i.e., 
Xn 7^ for all n. The corresponding orthogonal polynomials (ipn) are 



The second equation in (I4ip permits us to eliminate bn and formulate 
equivalently the recurrence only in terms of x„ and yn, 




(41) 



Ipn+liz) = {z + Xn)z'^ + Vn- 



1 2 I 
^n—l I Wn—l 



< 



^n—l 



(42) 



bn+l 



Vn-l 



Vn = 



^n—l 
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The second equation in (p2|) is solved by 



yn = (-!)" _ , (43) 



SO we only must care about the first equation in (j42p . 

If L = P + with P(z) = az + a G C*, /3 G M, we know that 

det^„(2;) = + (1 + — |yn.P)-z + a^n. oc = az^ + 2(3z + a. 

Therefore, 

Ol 1 + la^nl |2/n| (44) 

This implies that x„ = s^y^, G M, and the first equation of (|12|) is 
equivalent to 

x„ = 2(1; — — , uj = —. (45) 

That is, we have reduced the compatibility of (|4ip to the compatibility of 
(j45p for Xn, which can be rewritten in terms of s„ as 

s„ = 2u; ^^=^ (46) 

Sn-l |a| 

while the compatibility means simply that s„, 7^ for all n. If Sj 7^ for 
j < n but Sn = then the related solution is not quasi-definite but has only 
the first n + 1 MOP ^Ao, • • • , i'n- 

The key idea to calculate s„ is to write (|46p as a continued fraction 

s - 2a; ^ ' ' ^ ' ^ ' 

" 1 2a; 1 2a; |2a; \ sq 

According to the general theory of continued fractions (see for instance [27] ) , 

^ _ SpQn-l — Qn-2 
SO-fn.-l — Pn-2 ' 

where Pn and Qn satisfy the difference equations 

Qk = 2a;(5fc_i - Qk~2, Qo = 20;, Q_i = 1, 
Pk = 2a;Pfc_i - Pfc„2, i^o = 1, P-i = 0. 
Since Pi = 2a; = Qo, we get Qfe = Pk+i- 
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On the other hand, the recurrence and initial conditions for Pk show 
that Pfc = Uk{uj), where f/fc is the second kind Tchebyshev polynomial of 
degree k, 

^k+i _ x-{k+i) 

^i'('^) = 'x^x^i ' A = cj + Vc^2_i. 

The parameter A is one of the roots of the characteristic polynomial 

B{z) = A{-za/\a\) = z"^ -2UJZ + 1, (47) 

no matter which one because both of them are inverse of each other. 
Hence, 

Sof^n.-l(w) - C/„_2(w)' ~ 

As a consequence, the solution of the inverse problem is quasi-definite if and 
only if 

soUn{io)^Un-i{uj), n > 0. (48) 

In case sq Uj{oj) ^ Uj-i{oj) for j < n but sq Un{uj) = Un-i{uj), the related 
solution of the inverse problem is quasi-definite in P„ but not in Pn+i- 

Besides, from the solution for s„ we can obtain the rest of the variables 
of interest for the inverse problem. In particular, for n > 1, 



^n+l =yn = 



— \ n 



We can express these variables, as well as the quasi-definiteness condition 
(|48p . in terms of other parameters. For instance, following Algorithm 12, we 
can use as free parameters bi and xq. Then, using (|44p and the relation 

Uo = bi — xq, we get for some k G M*, 

a = KXo, /3 = I (1 - + 2Re(xo6i)) . (49) 

If we chose the approach of Algorithm 13, then the free parameters must 
be 6i and a,/3, so we should express sq, xq and yo in terms of them. Prom 
(|49]) . bearing in mind that k = \a\/ sq, we obtain 

_\a\ _a 1 A{-h) 

2 /3-Re(a6i)' ^° 2/3-Re(a6i)' ^° 2/3-Re(a6i)' ^ ^ 
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Finally, we can use the point of view of Algorithm II. This implies that 
we restrict our attention to the solutions of the inverse problem which are 
quasi-definite at least in an not only in Pi, which was the case till now. 
Then, according to Algorithm II, 61 and 62 could be used as free parameters 
too. This can be done using (jl^ and the relation 

yo xq- hi 
02 = yi = -— = — , 

Xq Xq 

which determines xq as the following function of hi and 62 > 

^0 = TrT9(^i +^1^2)- 

1 - \02\ 

Also, 62 can be expressed in terms of a, (3 and hi using (jSOp . which gives 

_ A{-hi) 
^ a(l-|6i|2)- 

The fact that the iterations (|46p generating the solutions of the inverse 
problem and the quasi-definiteness condition (|48p are given in terms of a, (3 
and So uniquely suggests the possibility of using these variables to param- 
eterize such solutions. However, this is not possible because an arbitrary 
value of a, (3 and sq can be associated with no value of hi or with infinitely 
many values of 61. Indeed, the first identity of (|50p can be written as 

1^1 - xq? = B{so), 
which shows that we have the following possibilities: 

• If B{sq) < there is no solution associated with a, (3 and sq. 

• If B{sq) = there is exactly one solution associated with a, (3 and sq: 
that one determined by a, (3 and hi = xq = SQa/\a\. 

• If B[sq) > there are infinitely many solutions associated with a, 
(3 and sq: those ones determined by a, /3 and any value of 61 in the 
circle with center xq and radius ^JB[sq). Therefore such solutions are 
parametrized by a phase. 

In consequence, given P{z) = az + (3, the inequality B{sq) > determines 
the permitted values of sq. The set of solutions associated with P and a 
permitted value sq will be called the circle of solutions for P and sq, and will 
be denoted C{P,so). Eventually B{so) = and C(P, sq) degenerates into a 
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single solution. The fact that the quasi-definiteness condition depends only 
on uj = P/\ct\ and sq means that all the solutions of C{P, sq) have the same 
number of MOP. 

It seems that the presence of the circles of solutions with similar prop- 
erties should have to do with some symmetry of the problem. The most 
obvious one is the rotation symmetry. If u = vL, then ug = vqLq for any 
angle 6, where the rotation of a Laurent polynomial / and a functional v 
are defined by fe{z) = /(e~*^z) and V0[f] = v[/_5i]. When ug = u we find 
that V £ Hi{u) implies vg G Hi{u). The only functional u which is invariant 
under any rotation is that one defined by the Lebesgue measure, so only in 
this case we can assure that Hi{u) is constituted by "circles of solutions" 
obtained by the rotation of one of them. 

Bearing in mind that we are identifying equivalent functionals and that 
the rotation of a functional preserves its quasi-definiteness properties, the 
rotation symmetry permits us to reduce the analysis of the set Hi (u) for the 
Lebesgue functional u to the case a = 1 because each "circle of solutions" 
has a representative with a monic polynomial P. However, the reduction of 
the analysis to such canonical cases is not possible for any other hermitian 
functional u. 

Nevertheless, the rotation symmetry of the Lebesgue measure is not re- 
sponsible of the circles of solutions C(P, sq) that we have found: the solutions 
of any circle C(P, sq) have a common polynomial P, while the solutions of 
a "circle of solutions" associated with the rotation symmetry are related to 
different polynomials P obtained by a rotation of one of them; furthermore, 
the rotation of a functional also rotates its Schur parameters around the 
origin, but the parameters bi of the solutions of a circle C{P,sq) are ob- 
tained rotating one of them around xq ^ 0. The search for the "symmetry 
transformations" relating the functionals of a circle C(P, sq) remains as an 
open problem. 

Some particular quasi-definite solutions deserve a special mention, i.e., 
the solutions with constant coefficients Xn, Vn, which are characterized by 
any of the statements of the following equivalence, which follows easily from 
the previous results, 

Sn = so,n > <^ Xn = Xq, Tl > <^ i/n = 0, n > <^ bn = 0, n > 2 <^ 
<^ 62 = <^ yo = <^ 6i = xo 44> A{-bi) = <^ A{-xo) = <^ B{so) = 0. 

Therefore, these constant solutions correspond exactly to the case where 
a circle of solutions degenerates into a single solution. The corresponding 
functionals are those ones associated with the Bernstein-Szego polynomials 
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^lJn+l{z) = {z + bi)z'^. Since —bi must be a root of A, such solutions can 
appear only when A has roots outside the unit circle, which corresponds to 
the Geronimus transformation of the Lebesgue measure. 

It is advisable to discuss the three possibilities (a), (b), (c) pointed out 
at the beginning of Section [3.2.21 according to the location of the roots of the 
polynomial A. The reason is that the qualitative behaviour of the solutions 
of the inverse problem depend strongly on the case at hand. Before doing 
this we must remark that, since B{z) = ^(— Aa/|a|), the roots Ci)C2 of A 
are related to the roots A, A~-^ of B through Ci = — Aa/|a|, (2 = —X~^a/\a\, 
and the three cases we want to discuss can be characterized in terms of to. 
Concerning this discussion, notice that, once 61 is fixed, any restriction on 
Lj becomes a restriction on the initial value xq by (|19|) . 

We will comment the asymptotics in each of the cases (a) , (b) , (c) using 
the notation p„ ~ qn to mean that lim(p„/(7„) = 1. 

{a) A{z)=a{z-C){z-l/C), C e ^ \io\ > I. 

This case corresponds to B having two different roots A, A~^ E M, 
thus we can suppose |A| < 1 so that ( = — Aa/|a|. Then, the quasi- 
definiteness condition (1481) becomes 



2n 



or equivalently 



1 — A 



1 _ |^|2n 

which can be also understood as a restriction on bi because, together 
with A, it determines xq through ([50]) . 

Given only a and /3, not any value of sq is permitted because B{so) 
can be negative. This happens when Ai < sq < -^2, where Ai, A2 are the 
roots A,A~^ of B but ordered so that Ai < A2. Therefore, the values 
of So associated with a solution of the inverse problem are those lying 
on (— cxD, Ai] U [A2,oo). Then, the corresponding sequence of MOP is 
infinite or finite depending on whether the quasi-definiteness condition 
(jSip is satisfied for any n or not. 

There are two quasi-definite constant solutions: = A, rE„ = — (" = 
bi, Un = and Sn = x„ = — 1/C = bi, 7/„ = 0. Both of them give 
rise to a Bernstein-Szego solution with 6„ = 0, n > 2, but the first one 
is positive definite with measure dm{z)/\z — ^p, while the second one 
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is indefinite. As we will see, the solution dm{z)/\z — is somewhat 
singular among the solutions of the inverse problem, so in what follows 
we will consider only sq 7^ A, i.e., xq 7^ —C,- Then, 

bn+l=yn^02 — C = -a=f -C )C 

xq + C abi + aC 

lim5„ = limy„ = 0, lims„ = A~^, limx„ = — 1/C- 

Furthermore, the related orthogonal polynomials obey the asymptotics 

^|;n+l{z) ~ -a^^^a - IC^C'S \z\ < |C|, 
abi + aQ 

^Pn+l{z)^{z-l/C)z^, \z\>\C\- 

We observe that the parameters of the indefinite Bernstein-Szego 
solution provide the asymptotics of the parameters for all the solu- 
tions except for dm{z)/\z — (^p. Also, the indefinite Bernstein-Szego 
polynomials (z — 1/C)z^ yield the large z asymptotics of the rest of 
MOP which solve the inverse problem, with the exception again of the 
positive definite Bernstein-Szego ones {z — 

(h) A{z) = a{z-C)^, (eJ ^ 1^1 = 1- 

This is equivalent to state that B has a double root A = u;G{— 1,1}, 
which is related to C by C = — Aa/|a|. No quasi-definite solution with 
constant Xn can appear now, thus sq / A and xq 7^ — C for any quasi- 
definite solution. The confluent form of the Tchebyshev polynomials 
Un{io) = {n + 1)X^ yields the quasi-definiteness condition 

Tl 

SO + >^—-r, n>0, (52) 
n + 1 

i.e., 

n 

xo r -C — — T' " > 0, 
n + 1 

where, once hi is chosen, xq is fixed by ([501) with (3 = X\a\. 

If we fix only a and /3, then sq can take any real value because now 
B is non-negative on M. 



35 



We have the relations 



soUn{uj) - Un-l{uj) ~ (so - A)nA", 



1- |6iP 



bn+i =yn^ h — = -a^= , 

hm6„ = hmy„ = 0, hms„, = A, hmx„ = — C, 
and the asymptotics of the corresponding orthogonal polynomials is 

ipn+iiz) , \z\ < 1, 

abi + aC n 

V'n+l(2) ~ (2-C)2", kl > 1- 

We see that in this case there is a so well defined asymptotics for 
any solution as in (a). However, contrary to |u;| > 1, the asymptotics 
of the frontier case \uj\ = 1 defines no quasi-definite solution of the 
inverse problem. 

(c) A{z) = a{z-Ci)iz-C2), C1/C2, CfcGT ^ \u;\ < I. 

Now B has two different roots A,A G T so that Ci = — Aa/|a| and 
(2 = —Xa/\a\. The quasi-definiteness condition (jl8|) reads as 

soImA"+^ / ImA", n > 0, 

that is, 

^o(cr'-c2+')7^c2"-cr, n>o, 

which again can be considered as a constraint on bi due to (jSOp . 

Concerning the possible choices of sq when fixing only a and /3, any 
real value of sq is possible since B is now positive on R. 

Analogously to case (b), sq 7^ A, A and xq 7^ — Ci, — C2 for any quasi- 
definite solution. Writing A = e*^, 6 ^ Zvr, and sq — A = |so — A|e*''', 

TT I \ TJ ( \ I ,| Sin((n + l)g + 7) 

sinC7 

thus the quasi-definiteness condition can be stated as 

n6' + 7^Z7r, n > 1, 
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and we find the identities 

sin((n + 1)6' + 7) ^ sin 6* 

= w /) , ^ = cos + - — , 

sin(nt^ + 7j tan(nC7 + 7) 

Q \ " j/o sin g 

o^l/ |so — -^1 sin(n0 + 7) ' 

which show that in this case Sn and \bn\ do not converge for any quasi- 
definite solution. 

The algorithm (j46p giving the solutions of the inverse problem for the 
Lebesgue measure can be interpreted as a Newton algorithm to find the 
zeros of a function. It is instructive to discuss the different behaviour of 
the associated Newton algorithm depending on the values of lu and sq. This 
approach sheds light on the different asymptotics found in cases (a), (b) and 
(c). Since we will discuss the behaviour depending on the values of u! and sq, 
we remember that, given P, there is a set of permitted values sq and each 
choice of sq determines a circle of solutions C{P, sq) which degenerates into 
a single solution when sq is a root of B. Remember also that the solutions 
of such a circle have the same number of MOP. 

The Newton algorithm for a real function f{s) of a real variable s is 
given by the iteration 

_ _ /(Sn-l) 

/ (Sn-l) 

Comparing this with (I46p we see that the algorithm providing the parameters 
Sn of the inverse problem for the Lebesgue measure can be understood as 
the Newton algorithm for a function f{s) satisfying 

/(s) s 

Solving the above equation we find three cases (Ai, A2 are the roots of B): 



(a) |w| > 1 ^ f{s) 

(b) 1^1 = 1 ^ f{s) 



(c) |w| < 1 ^ f{s) 




1 



s- Ail^i 
s — uj\ exp 



u> 



uj — s 



\/B(s) exp ( : arctan ( : | | , 
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Figure 1: (Case (a) - quasi-definite circle of solutions) First values of s„ for lj = |, 
Ai = ^, A2 = 2, cr„ = 2^^rp^3Y, ■^0 = 5 ^ This value of sq generates an 

infinite sequence (s„) such that s„ — s- A J monotonically for n > 2. Hence, the 
solutions of the associated circle C(P, Sq) are quasi-definite. 




-1.0 -0.5 0.0 0.5 1.0 1.5 



Figure 2: (Case (a) - non quasi-definite circle of solutions) Values of s„ for a; = |, 
Ai = i, A2 = 2, So = = The iterations stop at n = 2, thus the solutions of 
the circle C(P, sq) have only the MOP ■0o, fpi, "02- Since the set {cr„} is infinite, 
there exist non quasi-definite solutions with an arbitrary number of MOP. 
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-2-10 1 2 3 4 

Figure 3: (Case (b) - quasi-definite circle of solutions) First values of s„ for uj = 
Ai = A2 = 1, cr„ = -p^, So = f ^ {cti}- The situation is similar to Figure 1, but 
now Ai = A2. 
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Figure 4: (Case (b) - non quasi-definite circle of solutions) Values of s„ for ui = 
Ai = A2 = 1, So = = |- The situation is similar to Figure 2 but now Ai = A2 
and we have chosen sq so that the solutions of the circle C(P, sq) have four MOP. 
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-2-10 1 2 3 

Figure 5: (Case (c) - quasi-definite circle of solutions) First values of s„ for w = |, 

Xi.2 = I ± |«, cm = 5j^^^||pTT, So = -^1 i {(Jn) C Q. The solutions of 
the associated circle C{P, sq) are quasi-definite because sq generates an infinite 
sequence (s„) which oscillates indefinitely around the origin. 
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Figure 6: (Case (c) - non quasi-definite circle of solutions) Values of s„ for = |, 
Ai.2 = I ± |i, sq = ^ ^ W" '^^^ iterations stop at n = 4, thus the solutions of 
the related circle C(P, sq) have only five MOP. Like in Figures 2 and 4, the set {(Jn] 
is infinite (but, on the contrary, (un) is not monotone neither convergent) because 
X\ 2 are not roots of the unity, so there exist non quasi-definite solutions with an 
arbitrary number of MOP. 
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-2-10 1 2 3 

Figure 7: (Case (c) - quasi-definite circle of solutions - periodic case) Values of s„ 

for Ai,2 = e±*f, a„ = Jg| f^^^ , so = I i = {0, A l/^/2, ex.}. 

Like in Figure 5, the solutions of the associated circle C(P, sg) are quasi-definite 
but, on the contrary, the sequences (s„) and (tT„) are periodic with period 4 because 




-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

Figure 8: (Case (c) - non quasi-definite circle of solutions - periodic case) Values 
of s„ for uj = Ai,2 = e^*Tj sq = a2 = Like in Figure 6, the solutions of the 
circle C{P,sq) are non quasi-definite, although in this case there exist only three 
MOP. Indeed, contrary to Figure 6, there is no non quasi-definite solution with 
more than three MOP because (t„ takes only three finite values: (Tq = 0, CTi = 
and (72 = l/\/2- The picture, which can be understood also as the inverse Newton 
algorithm starting at the origin which yields ((t„), shows clearly that ~ oo 
because the corresponding tangent line becomes any of the two asymptotes. 
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The typical behaviour of the iterations in these three cases is shown in 
Figures 1 to 8, which represent the function /(s) as well as some of these 
iterations for different choices of u and sq. In any case the function f{s) 
is analytic in M \ {Ai, A2} and has a minimum at s = which can stop the 
iterations, giving rise to a circle C{P, sq) of non quasi-definite solutions but 
with a finite segment of MOP with the same length for all the circle. 

When > 1 the function /(s) diverges to 00 at s = A and vanishes at 
s = A~^, where A is the root with smallest module among Ai, A2. Indeed 
A~^ is also the absolute minimum and, despite the visual effect in Figure 
1 at A-\ / G C(^({A-i}) so /'(A~^) = 0. Excluding the case sq = A, the 
iterations, which must start at a point of (— c«, Ai]U[A2, 00), always converge 
to A~^ (corresponding to a circle of quasi-definite solutions) or they stop at 
the origin after a finite number of steps (corresponding to a circle of solutions 
with only a finite segment of MOP). 

If w = ±1, then lims^;^^ /(s) = 00 and lims_>;^± /(s) = lims^;^± f'{s) = 
0, where A = Ai = A2 = ±1, which plays again the role of an attractor where 
the iterations converge (circle of quasi-definite solutions) while they do not 
stop at the origin (circle of solutions with a finite segment of MOP). 

On the contrary, f{s) has no divergence neither zero when |ti;| < 1, and 
the origin is then the absolute minimum. In this case, as far as the iterations 
do not reach the origin (circle of solutions with a finite segment of MOP), 
they oscillate indefinitely around such a minimum (circle of quasi-definite 
solutions). 

In any case, for each value of u>, the values of sq associated to non quasi- 
definite solutions can be obtained by the inverse Newton algorithm starting 
at the origin, so they form a sequence ((t„) given by 

an = , o-Q = 0. (53) 

If So = CTn, then sj = cTn-j 7^ for j < n and s„ = 0, hence the solutions 
of the related circle C{P, sq) have only n + 1 MOP. When > 1, ((T„) is a 
monotone sequence with limit A, but if < 1 then ((T.„) is non convergent 
and oscillates around the origin. Eventually, (T„-i = 2uj and the above 
iterations stop. To understand this fact notice that (08]) shows that 

"""" Unico) 

if Un{uj) / 0, otherwise an has no meaning because no value of sq can satisfy 
soUn{uj) = Un-i{uj) when C/„(a;) = 0. The recurrence for Un implies that 
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cr„_i = 2uj iff Un{u)) = 0, so this is exactly tlie case wliere (T„ does not exist 
and, besides, an+i = = (Jq, hence the values of crj, j > n + 1, are simply 
a reiteration of the values for j = —1,0, . . . ,n — 1 if we define it_i = oo. 
Therefore, (j53p always works for n > —1 if we assume that an-i = 2a; gives 
(7„ = oo, which leads to an+i = and yields a periodic sequence (aj) in 
M U {oo} with period n + 1. 

Summarizing, if w is a zero of Un, which can hold only when \uj\ < 1, 
there is a finite number of non quasi-definite circles of solutions C(P, sq), 
those ones related to the initial values sq £ Wj}^Zi- Furthermore, if n is the 
smallest index such that Un{uj) = 0, the quantities fjj, j = 0, . . . , n — 1, are 
different from each other, hence there are exactly n — 1 non quasi-definite 
circles C{P, sq), and the length of the corresponding finite segments of MOP 
runs from 2 to n when sq = ai, . . . , (Jn-i- Therefore, there are no non quasi- 
definite solutions with more than n MOP. 

On the contrary, if Un{io) / for all n, then aj 7^ ak for j 7^ k, thus 
an infinite denumerable set of non quasi-definite circles C(P, sq) appear, 
which correspond to sq E {o'j}j^Q- In this case, given any ra G N, there is 
exactly one non quasi-definite circle of solutions with only n + 1 MOP, which 
corresponds to sq = an- 

As a final remark notice that Un{(jj) = means A^"^^ = 1, A 7^ ±1. 
Therefore, not only the sequence (cj), but also {Uj{u))) is in this case periodic 
with period n+1, so (sj) shows such a periodic behaviour too no matter 
the choice of sq. 

4 Applications of these techniques 

The characterization we have obtained for hermitian functionals related by 
polynomial perturbations is not only interesting by itself, but provides an 
efficient tool to answer different questions concerning orthogonal polynomi- 
als on the unit circle. In this section we will show two examples of this. The 
first one exploits the fact that a polynomial perturbation is equivalent to 
a linear relation with polynomial coefficients between two sequences of or- 
thogonal polynomials and their reversed ones. The second one deals with a 
problem concerning associated polynomials, which can be solved due to the 
formulation of a polynomial perturbation in terms of a difference equation 
for two sequences of Schur parameters. 
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4.1 Orthogonal polynomials and linear combinations with 
constant polynomial coefficients 

There are in the hterature different resuhs on the orthogonaUty properties 
of Unear combinations of orthogonal polynomials. In particular, it is known 
that, if {^n) and (V'n) are MOP on the unit circle, a relation like 

r 

Ipn+r = y^(Aj>y?n+i + ^j.nV'n+j), >^j,n, Hj,n G C, Ao,n / 0, 71 > 0, (54) 
j=0 

forces (^n) to be Bernstein-Szego polynomials when r > 1 (see j20j). The 
result is so strong that it holds assuming (j54p only when n > uq for some 
no, and even if we suppose that the sum in (|54p is up to and index r(n) 
depending on n, with the simple restriction 1 < r(n) < n/2 for n > uq (see 

EH). 

A way to escape from this triviality is to consider a more general relation 
than dM])- Identity §^ imphes that tpn+r G (^Pn-2)^"+'' C (z^Pn_,._i)^"+'- 
for r > 1, where the orthogonality is understood with respect to the func- 
tional associated with {(pn)- Thus, Lemma shows that is a particular 
case of 

tjjn+r = XniPn + Yn^p^, X„ G P^, y„ € ¥r-l, n > 0. (55) 

However, contrary to (j54p . a relation like (j55p can hold for non trivial MOP 
{(Pn) and (^n)) since it is always equivalent to a polynomial perturbation 
relation between the corresponding orthogonality functionals due to Theo- 
rem [2TT] and the subsequent comments, together with Proposition 12. 141 the 
hermitian functionals u and v associated with {(fn) and {ipn) must be related 
by ti = vL where L = P + is given by a polynomial P with degP < r; 
the condition X„(0) 7^ 0, which holds for no n or simultaneously for all n, 
characterizes the case degP = r. 

In this section we will show that the freedom enclosed in (j55p is large 
enough to yield non trivial solutions even when imposing very strong condi- 
tions on Xn and Y^. More precisely, we will find all the pairs of sequences of 
MOP (^n) and {'Pn) related by (155p with constant polynomials coefficients, 
i.e., 

^l^n+r=X^n+Yipl, X G P„ F G P^-i, n > 0. (56) 

This is not only an academic problem, but its importance relies on the 
fact that the constant solutions should play the role of fixed points with 
respect to the asymptotics of the polynomials X„, Yn related to the quasi- 
definite solutions of Hr{u). Therefore, some of these fixed points should act 
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as attractors whose study could give information about the asymptotics for 
the quasi-definite solutions of Hr{u), similarly to what happens in Example 
[3X21 

Relation (j56p can be rewritten, together with its reversed, as 

^n+r = ^$n, = ( ^y* ) ' ™ ^ 0' 

and the polynomial perturbation is recovered by j4 = detX. 

As follows from Theorem 12.131 and Proposition 12.141 the problem we 
want to solve is equivalent to the recurrence T^+^-A' = XSn, n > 1, and the 
initial condition X^q = i.e., 

anX - bn+rX* = {z - 1)Y, " ' (57) 



ipr = x + Y, XeFr, Ye 



r-1- 



r 



If y = 0, equations ([57|) yield bn+rX* = UnX and = X. Since V'; 
and ip* have no common roots, we find that a„ = 6n+r = for n > 1. 
This situation corresponds to u being the functional associated with the 
Lebesgue measure and MOP <^n{z) = z'^, and v a Bernstein-Szego type 
functional with the first r + 1 MOP generated by arbitrary Schur parameters 
bi, ... ,br £ C\T, while ipn+riz) = z'^tpr{z) for n > 1. 

Let us find now the solutions with Y ^ 0. Denote for convenience a = 
and b = bn+r- The first equation of ([57|) simply says that Y is proportional 
to a self-reciprocal polynomial in Pr-i and |6| = \a\. Using such equation 
and bearing in mind that V'r = X + Y and tjj* = X* -\- zY*, we can eliminate 
X and X* in the second equation of (I57p . which becomes 

aV'r - bip; = [z{l - a) - (1 - a)] Y. (58) 

Therefore, 



^^"^ - T^m) ^^c " ^^^^ 



1-a VV'r(C) 

where we have used the Christoffel-Darboux formula for the n-th kernel 

Kn{z,C) = E?=o^7Vi(^)V'j(C) associated with the MOP (ipj). 
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As a consequence, given tpr, the solutions of (j57p are determined by an 
arbitrary choice of a € C \ T: (j59p provides b and Y self-reciprocal in Pr-i 
up to a factor, solving the first equation of (j57p . and finally X = '(/',. — y 
solves the second equation of ([57]) . 

On the other hand, given X, Y, let us see how many solutions a, b of 
()57p we can expect. If we suppose two different solutions a, b and a', b', (j57p 
gives 

na-a')Y = {b-b')Y*, 

\{a-a')X = {b-b')X*. ^ ^ 

Then, Y* oc Y, X* oc X and, using again (j57p . we find that y = or 
X oc (z — l)y. In the first case ipr = X, which is not possible because 
X* X. In the second case Y divides ipr = X + Y, which implies that Y 
is a constant because Y* oc Y. Hence, X{z) = z — 1 and the polynomial 
modification must be of degree r = 1. 

As a conclusion, given X, Y, the equations (j57p have at most one solution 
a, b when the degree r of the modification is greater than 1, or when it is 
equal to 1 but X{z) ^ z — 1. Thus, concerning the MOP related by (j56p we 
have to distinguish two cases depending on the degree r of the modification. 



• r > 1 

In this case, given X, Y , the Schur parameters a^, ^^n+r must be 
constants of equal modulus for n > 1: the unique solution a, b of 
equation (j57p . Furthermore, for any choice of a, 6i, . . . , 6^ G C \ T the 
system (j57|) has a unique solution m. X,Y ,b obtained through (j59]) and 
the relation X = if^^ — Y . In other words, the MOP related by (j56p are 
those corresponding to a sequence of constant Schur parameters 
(a, a, . . .) and those (V'n) related to a sequence (6i, . . . , 6^, b,b, . . .) of 
Schur parameters, where a,bi, . . . ,br £ C \ T are arbitrary and b is 
given by (j59p . The MOP related by (|56p are thus parametrized by 
a, 6i, . . . , 6r S C \ T. 

• r = 1 

If 7^ z — 1 the conclusions are similar to those corresponding to 
r > 1. However, when X{z) = z — 1 the system (I57p has infinitely 
many solutions no matter the choice of y = y G C. To see this, let us 
write ([57l) explicitly, 

' ay = by, 
< a + b = y, 
^bi = y-l. 
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Since 61 ^ T forces y ^ 0, the solutions a, b are all the symmetric points 
of the perpendicular bisector n(y) of the segment [0,y]. Therefore, 
the solutions corresponding to X{z) = z — 1 can be construct in the 
following way: choose bi S C \ T, which determines y = 61 + 1; for 
each n > 1 choose a„ G n(?/) \ T and bn+i G n(y) as its symmetric 
point with respect to the segment [0,y]. This procedure generates all 
the sequences of Schur parameters (bn) whose MOP {^n), (V'n) 

are related by 

Hence, the solutions with X{z) = z — 1 are parametrized by 61 G C\T 
and an infinite sequence (ai, 02, . . .) lying on 11(1 + bi) \ T. 

On the other hand, the solutions with X{z) 7^ z—1 are parametrized 
by 61,0 G C \ T with a ^ 11(1 + 61), and the corresponding pair of 
sequences of Schur parameters is given by (a, a, . . .) and {bi,b,b, . . .) 
with b = a{C + bi)/{l + biC)- This yields ah the MOP related by 

ipn+i{z) = {z + x)ipn{z) +yip*n{z), X, y G C, X -1. 



Moreover, from this equality for n = and (I59p we find that the 
parameters x, y related to a choice of 61 and a are 

«(l-|&iP) 

x = bi-y, y = - - -. (61) 

(1 - a) + 61(1 - a) 

Concerning the possible values of the polynomials X and Y, we have 
to point out that Y must be proportional to a self-reciprocal polynomial in 
Pf—i, as follows from ()57p . Indeed, ()59p shows that Y{z) is proportional to a 
kernel Kr-i{z, for some C £ 1", thus it has exact degree r — 1 unless Y = 0. 
On the other hand, X is a monic polynomial of degree r which can not be 
proportional to a self-reciprocal one unless r = 1 and X{z) = z — 1, as follows 
from the reasoning in the paragraph after (I60p . This, together with the fact 
that ijjr = X+Y must be an orthogonal polynomial, are necessary conditions 
which must be fulfilled by the polynomial coefficients X, Y. Nevertheless, 
they are not sufficient conditions for the existence of MOP satisfying (j56p . 
To see this consider the case r = 1, where these conditions become 

X{z)=z + x, Y{z)=y, xGC\TU{-1}, x + yeC\T. (62) 

However, solving (j61l) for bi and a we get 

bi=x + y, a = y- r-rir, 

1 — xr 
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which shows that to get the alluded necessary and sufficient conditions for 
r = 1 we must add to (j62p the following one 



\y\ / -r- — if kl / 1- 

l + x 

Concerning the polynomial perturbation L = P + such that u = vL, 
we know that A oc XX* — zYY* . Hence, when X{z) = z — 1 we find 
that P{z) (X z + (|yP/2 — 1). As for the rest of solutions, related to Schur 
parameters (a, a, . . .), {bi, . . . ,br,b,b, . . .) with b given in ([59]) . we only know 
that deg-P < r. The inequality degP < r is characterized by any of the 
statements of the following equivalence, which follows from the previous 
results and the recurrence for (tpn), 

degP <r^ X{0) = 0^ Y{0) = br ^ b = br ^ 
'^"r = a <^ br = a— — -—. 

That is, among the values of a,6i, . . . ,6^ which parametrize the solutions 
with X(z) ^ z — 1, the inequality degP < r holds for those ones with br 
determined by o, 6i, . . . , br-i through br = aipr-iiO/i'r-iiO- The solutions 
with degP < r correspond to bn = b for n > r, while the solutions with 
degP = r are those ones with given by . . . ,br,b,b, . . .), br 7^ b. 
Notice that each solution with degP < r has a sequence with the form 
(bi, . . . ,bs,b,b, . . .), bg 7^ b, for some s < r, and then degP = s and one 
can find new polynomial coefficients X £ ¥s, Y £ Ps-i such that ipn+s = 
Xifn + yy^n^ n > 0. In any case, b = a'0j(C)/'0^(C) for 3 > degP. 



4.2 Associated polynomials and polynomial modifications 

Given a sequence (ipn) of MOP with Schur parameters the associated 
polynomials are those MOP (ipn) with Schur parameters «n = ^n+i- 

Despite the similarity of their Schur parameters, the corresponding orthog- 
onality functionals can be quite different. We will consider the following 
question concerning such functionals: when is the functional u of the associ- 
ated polynomials {(fn) a polynomial modification of the functional v related 
to the original MOP (ipn)'^ We will answer explicitly this question for a 
polynomial modification of degree 1. 

According to Theorem [27131 this is equivalent to the existence of matrices 
Cn G Ji such that C„fi„+i = AnCn-i, Bn+i = An, n > 1, with Co e 1^°^ 
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satisfying the initial condition Cq^i = A^q. Let us denote P{z) = az + (3, 
a G C*, /3 G M. The recurrence for C„ can be written as 

aZ + Cn dn \ { I an\ ( I an \ f az + Cn-l zdn^i 



zdn CnZ + a J \an 1 / 1 / V C„_iZ + a 

for some coefficients c„ E M*, (i„ G C. Sphtting this matrix recurrence gives 
the equivalent system of equations 

Cfi -\- CLfid-n C-n — l ~l~ '^n^n— 1) 

a„c„ + dn = aan, (63) 
aan = a„c„_i + 

Taking determinants in the matrix recurrence and setting z = 0, we find 
that c„ = c„_i for n > 1, so Cn = cq for n > 0. Therefore, ([63]) reads as 

Oindn — 0,ndn—l, 

< an {a- Co) = dn, (64) 
, a„ (a - Co) = dn-i, 

although the first equation is a consequence of the others. 

Assume that a = co- Then, dn = for all n and the initial condition is 
A = a{z + l)(z + 6i), which is not possible because A is self-reciprocal while 

7^ 1. Hence, a 7^ co and the solution of (I64p is 



a„+i = A"ai, dn = A"(a — co)ai, A = -, n> 0. 

a — Co 

Besides, the initial condition 

az'^ + 2Pz + a = {az + co) {z + bi) + do(biz + 1) 
yields the parameters of the polynomial perturbation, 

a = 61 Co + do, (3 = ^{abi + bido + cq) = y(l - |6ip) + Re{abi). 

Taking into account that do = {a — co)ai, we can express a, /3, A, do, in 
terms of oi, 61, cq, 

ai{bi - ai) + (61 - ai) 



a = cq- 



/? = Co 



/ 1(1 _ |fej2x ^ Re[(ai(6i -oi) + (6i -ai))6i] | 



_ (6i-l)+ai(6i-l) 
(6i-l)+ai(6i-l)' 

(6i-l)+ai(6i-l) 

flQ — CqCi j p5 

1 — air 
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The fact that deg P = 1 means that a 7^ 0. This only excludes the possibility 
ai = bl, which gives A = 1 and thus corresponds to the trivial case a„ = bn 
for all n, i.e., u = v. 

Therefore, the arbitrariness in cq £ M* is simply the freedom of the 
polynomial perturbation in a multiplicative real factor, and the solutions of 
the problem are parametrized by ai,6i € C \ T with oi 7^ 61: the MOP 
(ipn) whose associated ones come from a polynomial perturbation of 
degree 1 of the orthogonality functional of (ipn) are those ones with Schur 
parameters (61, ai, aiA, aiA^, . . .), where A E T is the square of the phase of 
(61 — l) + ai(6i — 1). The associated polynomials have Schur parameters 
(oi, aiA, aiA^, . . .), so they are obtained by a rotation (pn{z) = X^cpni^z) of 
the MOP ((/)„) with constant Schur parameters (ai,ai,ai, . . .). 

We can use a, (3 and bi as free parameters too. The initial condition can 
be expressed as 

1 bi\ f co\ ^ f P-abi 
bl 1 J [do] V a 



with solutions 



This gives 



d-Reiabi) ^ Ai-bi] 
Co = 2 — 



ai 



l&iP l-l&il 



do A{-bi 



a -Co a(l-|6i|2)-2(/3-Re(Q6i))' 



a(l - Ifoip) - 2(/3 - Re{abi)) 
~ a(l-|6i|2)-2(/3-Re(Q6i))' 

providing a solution whenever cq 7^ 0, a and |ai| 7^ 1, i.e., 

Ai-bi] 



(5 7^ Re(a6i), |(1 - \bi\'^) + Ke{ab 



a{l - |6i|2) -2(/3-Re(a6i)) 



/I- 
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